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Abstract 

V 

O i Herein we develop a dynamical foundation for fractional Brownian Motion. A clear re- 

O ■ lation is established between the asymptotic behaviour of the correlation function and 

> ' 
• I— I . 

^ . diffusion in a dynamical system. Then, assuming that scaling is applicable, we establish 
a connection between diffusion (either standard or anomalous) and the dynamical indica- 
tor known as the Hurst coefficient. We argue on the basis of numerical simulations that 
although we have been able to prove scaling only for "Gaussian" processes, our conclu- 
sions may well apply to a wider class of systems. On the other hand systems exist for 
which scaling might not hold, so we speculate on the possible consequence on the various 
relations derived in the paper on such systems. 
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1 Introduction 

A stochastic process is in general characterized by two quantities, the probabihty density 
describing the random nature of the fluctuations and the correlation function describing 
how a fluctuation at a given time influences subsequent fluctuations. The statistics and 
the spectrum of the random variations are therefore independent quantities and both are 
necessary to completely describe a stochastic process. For example, a Gaussian distribu- 
tion may specify the statistics of a process, but the time dependence of the variance, i.e. 
the width of the distribution, depends on the correlation function or spectrum of the pro- 
cess. Given that in many physical systems the Gaussian distribution is a straightforward 
consequence of the Central Limit Theorem, investigators are often satisfled with measure- 
ments of the correlation function to describe natural stochastic processes particularly in 
the observation of large scale phenomena such as in geophysics. One such scientist who 
was most successful in this regard was Hurst, who was particularly interested in problems 
of Hydrology and the Nile river ^ . 

Hurst |l|, H, 1^ developed a method called rescaled range analysis, in which the span of 
a random process is devided by its variance, resulting in a new variable that depends on 
the time length of the data record in a particularly interesting way. Let us define the time 
average of the random process ^{t) over the interval of time r: introducing t, the discrete 
integer-valued time at which the observations are recorded, and r, the total integer valued 
time-interval considered, we have 



(0. = -E^W 




t=i 
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Let us also define X(t), the "accumulated departure" of ^(t) from the mean {^)t, i-e. 

X{t,r)^j:au)-{Or (2) 

u=l 

SO that the span of the process is defined by 



Sir) = max X{t, r) - min X{t, r) (3) 

1<£<T l<i<r 



Finally, let us also consider the standard expression for the variance 



V{t) 
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(4) 



t=l 

The rescaled Hurst analysis consists in studying the properties of the ratio 

R{t) = S{t)/V{t) (5) 

He found that for the time records of over 850 phenomena R is well described by the 
following empirical relation 

R{t) = {T/2f" (6) 

Here we refer to Hh as the Hurst exponent (Hurst used the symbol K for this exponent). 

Mandelbrot and co-workers @, |], H modeled this heuristic result using the theory of 
fractional Brownian motion. They showed that fractional Brownian motion would provide 
an explicit statistical realization of ( and that the theory would imply a reasonable 
interpretation of the parameter Hh- We emphasize that this important result supports 
the interpretation of natural phenomena in terms of fractal functions. However, this 
interpretation does not take into account the fact that fractal functions, as important as 
they may be, are still idealizations of natural phenomena. These idealizations are not 
applicable to all time scales. In this context we mention the recent efforts to establish the 
physical foundation of classical statistical mechanics using the concept of chaos R ||, 0] . 



These efforts rely on there being a wide time scale separation between the microscopic 
and the macroscopic dynamical regimes. The microscopic quantities of motion are valid 
on a very short time scale, where the conventional statistical arguments are inapplicable. 
Then, upon increasing of the time scale considered, as a result of the action of chaos, 
the system of interest exhibits those statistical properties which are associated with the 
conventional ideas of a canonical distribution and Gaussian statistics. 

The main purpose of the present paper is to provide a dynamical basis for the Hurst 
rescaled range analysis. We show that the theory of Mandelbrot and co-workers |^, |^ 
focuses on the asymptotic limit of the dynamical processes considered herein. 

There are a number of benefits arising from this change in perspective. First of all, 
adding a "dynamical dimension" to the Hurst analysis, can be profitably used to quan- 
titatively illustrate the evolution of a deterministic system away from a regular toward 
a totally chaotic process, the latter being, for many purposes, virtually indistinguishable 
from a stochastic one. Thus, the Hurst exponent may provide an independent means of 
distinguishing stochastic from chaotic motion [|1^]. Following this idea, we are tempted 
to speculate that the distinction between these two kinds of physical processes may be 
merely a question of time scale. 

Secondly we shall see that such a dynamical analysis, i.e., putting the fractal geometry 
of Mandelbrot into a dynamical setting, has the beneficial effect of rigorously establish- 
ing a connection between the Hurst coefficient and the behaviour of the autocorrelation 
functions at long times. We link the coefficient to the existence of positive or negative 
tails for the autocorrelation function of the random variable ^, 

m = mm) iv 
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In molecular dynamics, dynamical processes characterized by a long-time regime with an 
inverse power-law correlation function 

C{t) = ±k/r for t ^ oo (8) 



are denoted as slow- decay processes since the pioneering work of Adler and Wainwright |TT[1 
where ( |^) was first identified, correlation functions have been the subject of an intense 
debate |T^, From ( ^ we note that the power law decay can have either a positive 
or a negative tail. We shall see that the processes that Mandelbrot denotes as persistent, 
Hh > 1/2 in ( ^, are connected by our theoretical analysis to a positive tail, whereas those 
termed by him to be antipersistent, Hh < 1/2 in ( |^) , are connected by our theoretical 
analysis to the existence of negative tails. 

We support the analytic arguments presented in Section 2 and 3 with computer cal- 
culations done for a substantial number of different dynamical systems. The numerical 
results presented in Section 4 support the suggested dynamical approach to the "geomet- 
rical" theory of Mandelbrot, but they also suggest that the asymptotic time regime itself 
might be explored with different mathematical arguments, valid also for non-Gaussian 
statistics. 



2 Dynamical theory 

Let us now focus our attention on the following equation of motion 

X = e (9) 

The Hurst coefficient ( ^ was defined in terms of a discrete time process so any dynamical 
representation such as ( ^ ought to be discrete as well. However, to connect the process 
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with the field of molecular dynamics, we adopt a continuous time representation. The 
formal time integration of ( ^ yields 

x{t) = f' ^{t')dt' + x{0) (10) 
Jo 

Thus the variable x{t) undergoes a kind of motion with erratic fiuctuations induced 
by ^(t). Later on we shall define more precisely the nature of such "disordered" motion. 
For the time being we limit ourselves to the conventional language of statistical mechan- 
ics. Thus, rather than focusing on single trajectories we shall study the mean values 
of quantities like x"'{t). We make the simplifying assumption that the erratic variable 
^(t) fiuctuates around zero. Thus {x{t)) = (x(0)) where the brackets denotes an average 
over an ensemble of realisations of ^(t), as well as the distribution of initial state of x(t). 
We are therefore obliged to study the mean value of x^(t). By averaging x^(t) over the 
fiuctuations and initial conditions, we obtain 

{x'it)) = fdt' f dt"{i{t')i{t")) + 2 fdt\i{t')xm + (a;2(0)) (11) 
Jo Jo Jo 

We make the assumption that the second term on the r.h.s. of ( [T^) vanishes, hypoth- 
esizing no correlation between the initial value of x and C,- This assumption certainly 
holds true when there exists a large time scale separation between the dynamics of the 
fast variable ^ and that of the slow variable x. However, this hypothesis must be used 
with caution in applications, which often refer to situations of slow decay ( |^), implying 
an exceptionally extended memory. We then adjust the limits of the time integrals and 
obtain 

{x'it)) = 2 fdt' f dt" {i{t')i{t")) + {x'm (12) 
Jo Jo 

Under the assumption that the process ^(t) is stationary, i.e., its moments are inde- 
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pendent of the time origin, so that 

{amt")) = {at'-t")m) (13) 

we derive from ( |12]) the following integro-differential equation of motion 

'x\t)) = 2 tmm)dT (14) 



dt Jo 

Clearly, the {x^{t)) appearing in the l.h.s. of ( p^ ) must be connected to the long-time 
diffusional regime described by 

(x2(i)) =Xi'^^. (15) 
It is evident that the physical bounds on the possible values of Hd are given by 

0<Hd<1; (16) 

= defines the case of localization, which is the lower limit of any diffusion process, 
and Hd = 1 obviously refers to the case of many uncorrelated deterministic trajectories, 
with x{t) — x{0) linearly proportional to time for each of them. The bound Hd < 1 is a 
consequence of the fact that a diffusional process cannot spread faster than a collection 
of deterministic trajectories! Finally, the condition Hd = 1/2 is obtained for simple 
Brownian motion, where the variance increases linearly with time. 

On the other hand, using the definition of the correlation function given in ( 0), ( p^) 
becomes 

2D{t) ^ -{x\t)) = 2 C{t')dt'. (17) 

We can now show that, using ( |l^, the deviation of Hd from the conventional diffusion 
prediction Hd = 0.5 can be explained if the correlation function C{t) exhibits a slow 
decay. The joint use of ( H), ( 0) and ( leads to the following long-time prediction 

— (x^(t)) = 2HDi2HD - l)Kt'"--' ~ 2C{t) = ±- (18) 



having assumed that the long-time hmit of the correlation function C{t) is dominated by 
the inverse power law of ( |^). The positive (negative) sign refers to the case of the solid 
(dashed) line in Fig. 1. 

From ( |TH|) we determine that Ho and a satisfy the relation 

Hd = 1- a/2 (19) 

obtained by matching the time dependences. It is also clear from the coefficients in ( |18|) 
that Hd > 1/2 implies a positive long-time correlation, whereas Hjj < 1/2 implies a 
negative long-time correlation. Let us now summarize the result of this simple theoretical 
analysis, with an eye to Fig. 1: 

Case exemplified by the solid line: 1 > Ho > 1/2; 1 > a > 
Case exemplified by the dashed line: < Ho <l/2; 2>a>l 

Thus, we see that the solid line correlation function of Fig. 1 leads to a superdiffusive 
behavior ranging from the standard diffusion {Ho = 1/2) to the ballistic behavior {Ho = 
1). The dashed line correlation function of the same figure leads to a sub diffusive behavior 
ranging from the standard diffusion to no motion at all. 

It must be stressed that for superdiffusive correlation functions ( [1^), D{oo) = oo, 
whereas in the case of subdiffusive correlation functions, D{oo) is finite or even vanishing. 
In this latter case we obtain what we define as classical Anderson localization. At early 
times in the diffusional process the mean square value of x{t) increases. Then when the 
negative part of the correlation function C{t) becomes important, the rate of diffusion 
decreases. When the negative tail completely compensates for the positive part of the 
relaxation process, the rate of diffusion virtually vanishes. At this late time stage further 
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diffusion is rigorously prevented and the diffusing particle becomes localized. Processes 



of this kind have recently been discovered |]14[ and the theory presented here affords a 



remarkably straightforward explanation of them. It is interesting that such processes 
should admit such a simple interpretation. 

We point out that in the simple theory presented in this section, the only significant 
assumptions made are the stationary property of ( |13|) and the absence of correlation 
between a;(0) and ^. No assumption was made on the nature of the statistics of the 
stochastic process x except that it has a finite correlation function. In the next section 
we shall rederive ( |lj) under the more restrictive assumption that the process ^ and 
therefore x is Gaussian. We also note that the standard case Hd = 0.5 is compatible with 
almost any kind of relaxation process. The only condition to fulfill is that the correlation 
function C{t) is square integrable over the time interval . Thus, if we exclude the case 
Hd = 0.5, we must invoke an inverse power law decay to explain the behavior given by 



( 15). This is so because ( |T^) implies the existence of "stationary" behavior and which, 
in turn, implies the existence of an inverse power-law decay. A power-law decay is the 
only way of "killing" the possibility itself of defining a time scale, and this, in turn, is 
a condition essential to explain the "stationary" nature of the diffusion regime of ( |^) 
with Ho 7^ 0.5. What do we mean by "stationary"? Upon increase of the time scale 
considered, the diffusion process is increasingly dependent on the tail of the correlation 
function, until it becomes totally dominated by the inverse power-law tail. In this regime, 
since the power-law decay of the correlation function implies no time scale is dominant, 
the diffusion process becomes stationary. The concept of a stationary diffusion process 
can be easily expressed by referring to the Hurst coefficient rather than to H^,. We shall 
see that Hh is time dependent, and that it usually reaches a stationary value for t ^ oo. 



This definition of " stationary" diffusion beliavior implies tliat tlie diffusion coefficient Hn 
can be identified with the asymptotic value of the Hurst coefficient H^- We shall see that 
this is frequently the case even if we shall only be able to rigorously prove it for Gaussian 
statistics. 



3 A Fokker-Planck treatment 

The next step is to determine the connection of the Hurst coefficient with the diffu- 
sion coefficient Hd- To do this we need to derive suitable expressions for the quantities 
appearing in ( for any given dynamical system. In part we follow the strategy of 
Mandelbrot: the idea is that for a white Gaussian process it is simple to carry out the 
theoretical analysis, hence it is only necessary to find, for a given dynamical system, the 
corresponding "Gaussian" approximation in the appropriate "reduced" time scale. 



The details of the approach are given in reference |T5[. For the given dynamical 
system we want to study, we replace the dynamical equations with the following equation 
of motion for the probability density 

^^p{x,^,T;t) = Lp{x,^,T;t). (20) 

Here F stands for the entire set of variables necessary to describe the time evolution of ^. 
The " Liouville-like" operator L is divided into two parts as follows, 

L = Lj + Lb (21) 



where the "interaction" part determined by ( ^ yields and the phase space operator 

L, = -5^. (22) 
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and Lb defines the time evolution of the distribution of the variables ^ and F, determined 
by the dynamics of the corresponding set of variables. It is not necessary to define the 
explicit form of the latter operator since it depends in detail on the specific problem 
studied. 

We now use the Zwanzig projection approach |T6|, which consists in integrating the 
total distribution p over the degrees of freedom that are not of interest to us 

a{x,t) = Jd^dTp{x,^,T) (23) 



We apply this projection approach following the perturbation prescriptions of [14|, assum- 



ing that Li in ( ^2|) is a "weak" perturbation. Using this basic assumption, after some 



algebra described in detail in [|15|, we arrive at the result 



where 



-a{x,t) = E{t)—a{x,t), (24) 



S(t) = /*C(r)c/r (25) 



If ( ^4]) was used to determine the time evolution of {x'^if)) i.e., multiply ( by 
x"^ and integrate over x, it would reduce to ( In this sense ( |1^ and ( |2^) are 

equivalent. Unfortunately, the two equations are not truely equivalent because ( |17[) is 
obtained without making any assumption regarding the statistics of x, whereas ( is 
really the result of a second-order perturbation treatment, equivalent to assuming that 
the statistical process x is Gaussian. 

If now we rescale the time as follows 

t* = t^^°, (26) 
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can be written as 



Aa(x;r) = $(r)^a(x;r) (27) 



where, of course, 



, 1 I 

1 1 rt '^^D ijf'* I 

<^> in = I ^^^^ 

After some algebra it can be shown that 

hm = K (29) 

where is a finite constant. It is evident that in the asymptotic time hmit ( |2^) becomes 
a standard (time independent) Fokker-Planck equation, and the statistical process defined 
in terms of the scaled variable 

y = -^=. (30) 

becomes a Gaussian process with the distribution 

Piv) = ^= I -TT I • (31) 



2vr 

This result means that the random function x{t) is proportional to t*^^'^ or, in the original 
time scale, to t^°. Following Mandelbrot 0], we are then led to identify with Hh- 

This is the central theoretical result of this paper. In the asymptotic limit the dynam- 
ical process described by ( ^) becomes the fractional Brownian motion of Mandelbrot. 
However, in the short time regime the process can be substantially different from such 
a stochastic process. This is clearly illustrated by the numerical results of the following 
section. 

We note that Hd is the coefficient appearing in ( [T^). However, to be able to identify 
this coefficient with Hh it is necessary that we go through the Fokker-Planck equation 
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of ( ^4]), and this implies making a second-order approximation, absent in the derivation 
of ( |15|). In the next section we show, using numerical methods that Hh is very close in 
value to Hj;) even in cases when ( should not be valid. This suggests the possibility 
of establishing the property Hh = Hd without using the second-order approximation 
necessary for the derivation of ( p^). 

4 Numerical simulations 

We now present numerical simulations of a number of different dynamical systems, to 
which the Hurst analysis is then applied. Comments on the correspondence of Hd and 
Hh are discussed alongside the numerical data. For clarity, the different models are 
presented within separate subsections. 

4.1 Microscopic dynamics described by the Langevin equation 

The easiest way to create a stochastic process ^ with an exponential correlation function 
is through the Langevin equation 

i = ~im + f{t) (32) 

where f{t) is Gaussian white noise. Clearly, the relevant time scale here is given by I/7, 
and we expect that for times larger than I/7, we should find Hh = 0.5. This is confirmed 
(see Fig. 2) by the numerical simulations: the expected asymptotic regime is reached at 
shorter times for the case of larger 7. 

A similar model is given by the multidimensional Langevin equation 

i = wit) 

w = -Tw{t) -n^^{t) + F{t) (33) 
13 



In the case F ri, ( |3^) becomes indistinguishable from ( ^2|) with 7 = Vt^ /V. Again, as 
soon as the time scale considered is larger than the typical time scale of the system, we 
should have Hh = 0.5. However, here a problem arises. It is possible to consider either 
the time under the correlation function or the time over which the "quasi energy" in ( |33|) 
looses correlation: these two times are in general very different. By inspection of Fig. 3, 
it is clear that numerical simulations done for different Vt^ and the same 7 show a similar 
asymptotic behavior more or less in the same time region (compare the solid line and the 
dotdashed line in Fig. 3): whereas, when a different 7 value is considered (compare the 
solid and dashed lines), the asymptotic regime sets in at different times. This seems to 
imply that the relevant time scale is the time over which the energy looses correlation. 

Due to the linearity of ( ^3|), the Gaussian statistics of the stochastic force is transmit- 
ted to ^ and hence to x. Thus, the Central Limit Theorem is fulfilled and the asymptotic 
behavior must be characterized by Hh = Hd = 0.5. The parallelism between R and 
(a;^(t))^/^ of Fig. 4 shows that the transient behavior corresponds exactly to the time it 
takes for the second moment to reach the stationary condition, corresponding to standard 
diffusion. In the case of ( [33| ) this transient time is a complex function of the parameters 
F and Q, which we do not discuss here. But of course standard diffusion sets in over time 
scales which are connected to the time taken by the quasi-energy to decorrelate, hence 
justifying the behaviour shown in Fig. 3. 
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4.2 Bistable stochastic motion 

The next system we consider involves the motion of a stochastic particle in a bistable 
potential. The dynamics can be cast in the dimensionless form 

i = ^-e+f{t) (34) 

where f{t) is a stochastic Gaussian white force [see ( ^2])] of intensity D. The Hurst 
rescaled range analysis is applied to the variable ^. In the limit of small noise intensities 
D the system ( |3^) is characterized by two different time scales: one is the time scale 
of the inter- well relaxation time (Tr); the other one is the time scale of the intra- well 
dynamics {Tk, related to the Kramers rate for the system). For r ^ Tk the system 
should behave as a dichotomous random process, hence Hh = 0.5. Also, if applicable, in 
the range Tr ^ r ^ we expect Hh = 0.5: this is because before the Hurst analysis 
is able to single out the dichotomous random process (the behaviour for r ^ Tk ), the 
condition Tr<^t implies that, for the relevant r values, the dynamics is similar to that of 
( ^). This is confirmed by Fig. 5, where we show the result for R{t) obtained by digitally 
simulating ( |3^). For the simulations we use D = 0.1, which yields a Tk of approximately 
30, and a Tr, determined by the coefficients appearing in the force in ( of order 1. It 
is possible to note that for larger r values Hh approaches the correct value 0.5, and that 
for Th < T < Tk, Hh is smaller than it is, say, in the region r ~ Tk- Unfortunately, the 
condition T^i <^ t Tk is only weakly satisfied by our choice of the parameter D and it 
is not possible to observe exactly Hh = 0.5 in the range of intermediate r. 
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4.3 Noisy Lorenz model 

Another model to which we have apphed the Hurst rescaled range analysis is the Lorentz 
model perturbed by Gaussian white noise The perturbed Lorenz model is described 



by the stochastic differential equations 

i = ^i^-y) + fit) 

y = r^~y~^z (35) 
z = -z + ^y 

where f{t) is a Gaussian white process with intensity D. We set a = 8/3, 6 = 10 and 
r = 126.5, a set of parameters for which the Lorenz model is known to be periodic, as long 
as = 0. The chosen value of the control parameter r = 126.5 is roughly in the middle of 
a small periodic island within an r region for which the Lorenz model displays chaos. The 
addition of a small stochastic force (i.e., D ^ 0) "kicks" the system out of the periodic 
region, and leads to the observation of chaos [|l^. The transition periodic/chaotic motion 
as function of the noise intensity D is very smooth, and chaos, defined by a positive 
Lyapunov exponent, is observed for D > 10^^ — 10^"^. 

Applying the Hurst rescaled range analysis to the variable ^ of the Lorenz system, we 
expect that Hh will vary from around zero (the value for periodic motion) to 0.5 (diffusive 
dynamics) for increasing noise. Also, the r values for which the asymptotic behaviour is 
observed should become smaller as D is increased. The result of numerical simulations 
of ( ^) is shown in Fig. 6. Note that in Fig. 6 we have plotted -R(r)/r^/^ versus r, to 
more clearly show the asymptotic behaviour. For all curves, R{t) increases for small r, 
up to r ~ 1, which is the period of the noiseless Lorenz system for the chosen parameters. 
In the case of very small noise intensities (full and dotted line), we have that for larger 
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r, -R(t) goes like r°; behaviour typical of periodic motion. When we increase the noise 
intensity (small dashes line), after an initial decrease, R{t) increases: it eventually leads 
to Hh = 0.5, for still larger r values, but on the r range shown the asymptotic behaviour 
is not yet established. For this D value, we remark, that the Lorenz model is only very 
weakly chaotic. Finally, for even larger noise intensities (large dashes line) the departure 
of -R(t) from the noiseless curve takes place at yet smaller r values, and the expected 
asymptotic behaviour {Hh = 0.5, on the figure, represented by an horizontal line) for 
R{t) is clearly identified. 

4.4 A two-dimensional potential 

We next study the deterministic motion of a particle in the two-dimensional potential 

V{^, y) = cos U + yVs) + cos U - yVs) + COS (2y/V3) (36) 



which defines an infinite lattice of triangular symmetry (see ||18| for more details). The 
equations of motion we have integrated have the form 

i = V 

dV{^,y) 
^ = 7^ — 

y = w (37) 

^ ^ dVi^,y) 
dy 

It is known that when the energy of the system is efficiently small the particle moves 
from the bottom of each triangular cell, only occasionally wandering from cell to cell. As 
a consequence, the self-diffusion coefficient shows a peculiar behaviour as the energy is 
changed. 
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It is clear that in principle this motion could lead to anomalous diffusion: in particular, 
as the energy is decreased the "periodic" motion within each cell should become more and 
more dominant in the dynamics of the particle. Obviously, the "random" diffusion from 
cell to cell is still in place, hence we only expect some weak departure from Hh = 0.5 as 
the mechanical energy is decreased. By inspection (see Fig. 2 in |]18|) it is clear that the 
autocorrelation function of the velocity v has a visible negative tail (which is not at all 
surprising, remembering the importance of the periodic motion within each cell). We then 
expect that as the energy is decreased Hh should take on values smaller than 0.5. This 
is confirmed by the Hurst rescaled range analysis applied to the variable ^, and shown in 
Fig. 7. 

As clear from the discussion in the previous section, the fundamental question is 
whether the coefficient Hh should be related to if/) in situations of anomalous diffusion. 
We compare the quantity R{t) (Hurst analysis) and (f^(r))^/^ (see Section 2) in Fig. 8: 
we used an energy value equal to -0.90, for which the dynamics is supposedly anomalous 
{Hh = 0.43). The clear parallelism between the two curves at large times establishes that 
the diffusion is indeed anomalous as suggested by the Hurst analysis. 



4.5 Standard map 

We now present the results obtained in a discrete model, i.e. for the standard map 

Xt+i = Xt + — sm Ot 

ZTT 

Ot+i = Ot + xt+i (mod 27r) (38) 
The standard map is very convenient to test our interpretation of the Hurst rescaled 



analysis analysis: it has been recently shown ||T9[ that for appropriate parameter values 
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within the chaotic regime anomalous diffusion should arise. 



According to |T^, the anomalous diffusion is caused by chaotic orbits sticking to 
critical tori encircling accelerator mode islands. For this reason the correlation function 
C(t) should have the power law dependency of ( |18|). We studied the map for the same 
values considered in [|T^, i.e., for K = 3.86, 6.4717, 6.9115 and 10.053: the dynamical 
variable we considered for the analysis is the quantity = Xt+i —Xt. A standard diffusion 
behaviour is expected for K = 3.86 and 10.053, and an anomalous diffusion behaviour for 
the other two K values. That this is qualitatively the case is clearly shown in Fig. 9, where 
we have plotted -R(t) normalized to r^/^ versus r for the different fC's: as expected, the 
curve is an horizontal line (standard diffusion) for K = 3.86 and 10.053. We would now 
like to understand whether for the anomalous diffusion case we have some correspondence 



between the Hurst rescaled range analysis and the numerical work of |19|. In the case 



K = 6.9115, in |19| it is reported that theory and numerical simulations lead to 



« ~ ^ (39) 



Let us insert ( ^9]) into ( [T^). Adopting the notation of |T^, i.e., ( = 2Hd, we obtain 

( ~ 1 = 1.3333333... (40) 

According to [|1^ this prediction fits very well the result of the numerical calculation on 
diffusion. 

We remark once more that for K = 3.86 and 10.053, values for which the authors of [|T^ 



observe a standard diffusion we obtain, with a very high degree of accuracy Hh = 0.5; and 
where {K = 6.9115 and 6.4717) anomalous diffusion is predicted we obtain a value of Hh 
significantly different from Hh = 0.5. Moreover, for K = 6.9115 we obtain C = 1.2330, to 
be compared to C = 1.3333 from ( ^0]) . 
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Table 1 summarizes the situation. Note that the fourth column, (sim denotes the result 
of the numerical simulation (the only value there reported corresponds to K = 6.9115) 
and that the last column, ({Hh) reports the values of ( corresponding to the value of Hh 
evaluated numerically (third column). 

5 Conclusions 

It must be pointed out that the Mandelbrot analysis, leading to = Hd, is essentially 
based on the Central Limit Theorem, assumed to be valid even when anomalous diffusion 
occurs. The adoption of the Fokker-Planck treatment of Section 3 leads to the following 
time dependence of the x-distribution: 



It is then evident that the moment {x"'{t)) rescales in time as t"^^ . Since the Hurst rescaled 
range analysis refers to a quantity with the same dimension of x it is evident that it leads 
to 



The functional form ( |4T| ) suggests that in general, after an initial transient, the 
probability distribution a{x; t) should perhaps be described by the equation 



If the rescaling of ( applied, we would have that indeed in general Hh = H]j. 

Let us now briefly discuss some possible forms of F. There are three possible condi- 
tions: 




(41) 




(42) 




(43) 
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(i) /3 = 1/2, F is a Gaussian function of its argument. This is the standard diffusion 

process. 

(ii) (3 7^ 1/2, F is a Gaussian function of its argument. This is the fractional Brownian 

motion process. 

(iii) (3 7^ 1/2, F is not a Gaussian function of its argument. Note that this occurs for a 
Levy stable process pO |. 



It must be pointed out, however, that from a physical point of view it is hard to 
imagine a diffusion process with a deterministic origin agreeing with ( ^), and thus 
falling under case (ii), even in the case H > 0.5. The reason is that as established by the 
theoretical analysis of Section 2, this anomalous behavior comes from an anomalously slow 
correlation function, namely the correlation function of ( |^) with a < 1. In this physical 
situation, there is no hope to realize ( ^T]) as an effect of the Central Limit Theorem: 
the original process must be already Gaussian! In other words, if there existed Gaussian 
statistical processes leading to the slow decay of C{t), then the anomalous diffusion would 
be compatible with the time rescaling of ( ^). In our opinion, this is the physical nature 
of the fractional Brownian motion of Mandelbrot. It is the long-time asymptotic limit of 
a Gaussian process with an anomalously slow correlation function. 

We think that this situation might occur in statistical mechanics when the source of 
the Brownian motion, the statistical process x, refers to a physical condition characterized 
by a large number of degrees of freedom. However, in the last few years, there have been 
attempts to build statistical mechanics on chaos, without the joint action of a very large 
number of degrees of freedom 10) ^j; 01- In this physical situation x is a non-Gaussian 
statistical process and the Gaussian nature of diffusion stems from the action of the 
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Central Limit Theorem. If the process is not Gaussian, but it is fully chaotic, then the 
correlation function is exponential or, more generally, characterized by a well-defined time 
scale. 

In the special case when chaos and ordered motion coexist, however, the dynamical 
behavior of the system becomes extremely more complex, and a correlation function 
with an inverse power law might occur. This implies the breakdown of the time scale 
separation between diffusion and microscopic dynamics, and the consequent breakdown 
of the Central Limit Theorem itself. In this physical condition ( ^T]) cannot apply. 

Are there processes rescaling according to (iii), without implying the Gaussian as- 
sumption? We think that if an anomalous diffusion exists, then it is quite probable that 
it belongs to the class (iii). We are convinced that some of the processes examined in 
Section 4 belong to class (iii). If the rescaling in ( ^3]) with (3^1/2 holds true, then we 
conclude immediately that Hh = (3. However, this special condition raises the intriguing 
question of whether or not Hh = Hd, in this case. The dynamical realization of the dif- 
fusion process is expressed by ( Let us assume the X is the maximum possible value 
of ^. It is then evident that at the time N the ^ distribution must be contained between 
Xm and —Xm, with Xm = NX. Now, let us imagine that there are theoretical reasons to 
expect that the x distribution is characterized by long tails with an inverse power law 
1/x"^. It is then evident that the rescaling of ( ^3|) cannot apply to the whole space. This 
might generate a discrepancy between and Hd. Let us assume, for simphcity, that 
Xm = At. In such a case we get a rescaling of the same kind as ( ^) only for \xm/t\ < A. 
Thus the moments of the distribution rescale with a power law different from that leading 
to the time rescaling of ( ^31) . 

We wonder if a possible discrepancy between the two coefficients might be derived using 
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the data already available for the standard map [|T^. According to [|T9[, the distribution 
rescales as in ( ^3]), with /? = 3/5 for K = 6.9115. However, this distribution is truncated 
at the value |y| = |x/)f:^| = 1. If we assume that Hh is determined by the rescaling of 
( we obtain 

2Hh = 213 = 6/5 = 1.2 (44) 

thereby suggesting that the discrepancy between the numerical value 2Hh = 1.23330, 
obtained in this paper, and the numerical value 2Ho = 1.3333333, determined numerically 
by the authors of |T^, might be due not to the inaccuracy of the direct calculation of 



Hr) in (notice that the calculation of Hh is expected to more accurate than that of 
Hd), but it might rather depend on the breakdown of the condition Hh = H^, due to 
the non-Gaussian character of the distribution F of ( |43| ). 

We shall address these questions in further investigations. For the time being we 
must limit ourselves to saying that the Hurst rescaled range analysis seems to be an 
efficient numerical technique to explore how a dynamical system approaches its long-time 
asymptotic limit, or, equivalently, which is the short time dynamics of that asymptotic 
idealization referred to by Mandelbrot ^, ^ as fractional Brownian motion. 
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Table caption 

Table summarising the comparison between Hh from numerical simulations (third col- 
umn) as function of the parameter K (first column). The power of the correlation function 
tail, from the simulations, is shown under the heading a. For the definition of Q see text. 



Note that Qsim is the value computed in [0, and that C{Hh) is computed from our Hh 
values. 
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^ure captions 

Figure 1: Typical slow decaying correlation functions with a positive tail (solid line) 
and a negative tail (dashed line). 



Figure 2: the Hurst rescaled range analysis applied to the Brownian motion of q2. 
The labels indicate the resulting Hurst function R{t), for different 7. The straight 
lines are best fits to the appropriate power law. 



Figure 3: The Hurst rescaled range analysis is applied to the Brownian motion of 
The parameters values are: full line, = 10, 7 = 1; dashed line = 1, 7 = 10; 
and dot-dashed line = 1, 7 = 1. 

Figure 4: Comparison between the Hurst rescaled analysis (dashed line) and diffu- 
sion of the variable x (= (x^(r))-'^/^, see Section 2, solid line) in the system described 
by Parameter values chosen are = 1, 7 = 10. Note the parallelism between 
the different curves at large r times. 

Figure 5: The Hurst rescaled range analysis applied to the system described by 0: 
for D ~ 0.1 chosen, Tk ~ 30, and Tr ~ 1 (see text). The boxed numbers are the 
best fit Hh values evaluated in the region around the arrow head. 



Figure 6: The Hurst rescaled range analysis is applied to the Lorenz model ( pS]) . 
The curves are drawn for increasing noise intensities (solid line, D = 0, dotted line, 
barely visible in the bottom left corner,/} = 10~^, small dashes line, D = 10^^, large 
dashes line, D = 10~^). The quantity plotted is -R(r) normalized to r^/^ versus r, 
and in case of Hh = 0.5 we should have an horizontal line. 
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Figure 7: The Hurst rescaled range analysis is applied to the two-dimensional model 



of |36| and |3^. From top to bottom the energy (in round brackets Hu) has the value 
-0.60 (0.49), -0.80 (0.50), -0.85 (0.49), -0.90 (0.43) and -0.95 (0.39). We plotted 
the quantity -R(t) normalized to r^/^ versus r, and in case of Hh = 0.5 we should 
have an horizontal line. 

Figure 8: Comparison between the Hurst rescaled range analysis (solid line) and 
diffusion (= (t>^(r))^/^, see Section 2, dashed line) for the model of PU| and and 
for a mechanical energy equal to -0.90 {Hh = 0.43). The parallelism between the 
curve proves the relevance of the Hurst analysis in the calculation of the diffusion 
at large times. 



Figure 9: The Hurst rescaled range is applied analysis of the standard map, 
The quantity R{t) normalized to r^/^ versus r is plotted, and in case of Hu = 0.5 
horizontal line results. The different curves refer to different values of the parameter 
K (see text). We have: solid line, K = 3.86; small dashes line, K = 6.4714; large 
dashes line, K = 6.9115; and dot dashed line, K = 10.053. Standard diffusion is 
expected for K = 3.86 and K = 10.053. 
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Table I 



K 


a 


Hh 


Csim 


C{Hh) 


10.053 


0.9483 ±0.0058 


0.50317 ±0.0004 




1.0624 


6.9115 


0.80188 ± 0.00888 


0.61665 ± 0.00081 


1.3333 


1.2333 


6.4717 


0.88778 ± 0.00854 


0.57652 ± 0.00066 




1.15304 


3.86 


0.86957 ±0.00443 


0.50134 ±0.000307 




1.0268 
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